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ABSTRACT 

Globular cluster stars evolving off the main sequence are known to lose mass, 
and it is expected that some of the lost material should remain within the clus- 
ter as an intracluster medium (ICM). Most attempts to detect such an ICM 
have been unsuccessful. The Multiband Imaging Photometer for Spitzer on the 
Spitzer Space Telescope was used to observe eight Galactic globular clusters in an 
attempt to detect the thermal emission from ICM dust. Most clusters do not have 
significant detections at 70 /im; one cluster, NGC 6341, has tentative evidence 
for the presence of dust, but 90 /im observations do not confirm the detection. 
Individual 70 fim point sources which appear in several of the cluster images are 
likely to be background galaxies. The inferred dust mass and upper limits are 
< 4 X 10~^Mq, well below expectations for cluster dust production from mass 
loss in red and asymptotic giant branch stars. This implies that either globular 
cluster dust production is less efficient, or that ICM removal or dust destruction 
is more efficient, than previously believed. We explore several possibilities for 
ICM removal and conclude that present data do not yet permit us to distinguish 
between them. 

Subject headings: globular clusters: general — globular clusters: individual 
(NGC 104, NGC 362, NGC 1851, NGC 5272, NGC 5904, NGC 6205, NGC 6341, 
NGC 6752) — infrared: stars — stars: mass loss 
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Introduction 



Post-main-sequence mass loss affects botli present and future generations of stars. It 
determines tlie remnant state of stars of intermediate mass and below, as well as affecting the 
state and contents of the interstellar medium out of which new stars are formed. However, 
the rate and amount of mass lost by red giant branch (RGB) and asymptotic giant br anch 
(AGB) s t ars is a major uncertainty in stellar evolutionary theory. As lOriglia et al.l (120071 ) and 
CatelanI (120001 ) point out, existing formulations rely on empirical laws based on observations 
of Population I stars. There is little theoretical guidance on how mass loss should be included 
in stellar evolutionary models, or on how it affects population synthesis models based on the 
stellar models. 

Globular clusters (GCs) have long been the prime testing grounds for stellar evolution 
theory. With large populations of low-mass stars of the same age, metallicity, and distance, 
they should be ideal places to test ideas about mass loss. There are several observational 
confirmations of the theoretical prediction that stars more massive than 0.8Mq lose a sig- 
nificant fraction of their mass after leaving the rnain se quence. These include initial-final 
mass estimates of white dwarfs (^e.g. iKalirai et al.l |2008| ) and evidence from line profiles of 



red giants fe.g. jMeszaros et al.ll2008l ). To explain the horizontal branch it was argued (e.g.. 



Roodlll973l ) that much of this mass is shed on the RGB, but as AGE stars , post-AGB stars 



and PNe are also observed in GCs (lAlves et al.ll2000l : Ivan Loon et al.ll2007l ) some mass must 



stars which are currently losing mass, at rates of 10 ^ — 1 



20011 : lOridia et all 120021 : iBover et all 120061 : lOridia et all 120071 : llta et al. 



o yr M 


Ramdani & Jorissen 


Ita et al. 


20071). From these 



observations it is difficult to obtain an overall picture of time-averaged mass loss rates and 
duty cycles and their relation to cluster properties such as metallicity, since only a few stars 
at a time lose mass. 

A more global picture of mass loss within GCs might be provided by the intra-cluster 
medium (ICM). Mass lost from GC stars is expected to accumulate in the ICM until the 
cluster passes through the Galactic disk, when the ICM will be stripped by ram pressure. 
The ICM is expected to be mostly gaseous (although the ionization state is unknown), with a 
small amount of dust formed in the extended atmo spheres of red giants. The dust should b e 
optically thin and heated by starhght to 50-80 K (lAngeletti et al.lll982l : IPorte et all 120021 ): 
its thermal emission should peak in the mid-to-far-infrared (IR) and could be detected as 
an IR excess at these wavelengths. Although GC stars are clearly losing mass, detections of 
the ICM have been surprisingly difficult. Searches for intra-cluster gas have been conducted 



for several decades ( IBirkin shaw et al.lll983l : [Roberts 



1988 



Faulkner et al 



1991 



2001 : van Loon et al.i,2006i ): searches for dust ( Lynch fc Rossano 199ol : Knapp et al 



Freire et al. 



1995 
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Origlia et al.lll996l : iPenny et al.lll997l : iHopwood et al.lll999l : iMatsunaga et al.ll2008l ) are more 



Only one ( 


Glalactic GC,NGC 


Evans et al. 


2003 




Bover et al. 



20061 ). Diffuse infrared emission near NGC 6402 may be associated with the cluster or with 
Gala ctic cirrus, while NGC 5024 has an infrared point source which may be intracluster 
dust (IMatsunaga et al.ll2008l ). An a dditional cluster, NGC 6356, has a tentative detection of 
dust at submillimeter wavelengths (IHopwood et al.l 119981 ). All of the detections imply dust 
masses well below the expected levels. 



This paper presents results of a Spitzer (IWerner et al.ll2004j : iGehrz et al.l 120071 ) survey 
of eight Milky Way GCs usi ng 24 and 70 nm o bservations with the Multiband Imaging Pho- 
tometer for Spitzer (MIPS; iRieke et al.ll2004l ). We describe our observational strategy and 
data analysis techniques in §2. Section 3 introduces our results, including the detection of 
excess IR emission in several of the target clusters and upper limits on the remaining clus- 
ters. Section 4 explores the implication of these new mid-IR results, including a comparison 
with predicted dust masses. We conclude in §5 with a summary and suggested follow-up 
observations. 



2. Observations and data reduction 



The target list (Table [T]) comprises the Milky Way globular clusters that we believed 
to be the best candidates for ICM dust detection. These clusters are reasonably nearby, at 
high galactic latitude (to reduce confusing fore- and back-grounds), and have high escape 
velocities. They are among the most massive Galactic GCs. More massive clusters are more 
likely to harbor a larger number of stars in the short-lived superwind phase: if the ICM is 
being removed from the cluster on short time scales, then we are more likely to see it in a 
more massive cluster where the ICM material is being replenished more quickly. Clusters 
with escape velocities that are much larger than AGB wind speeds (~ 10 km s~^) will also be 
more likely to retain ICM material. At the time this program was planned, all ha d existing or 



scheduled Spitzer mid-IR observations with the InfraRed Array Camera (IRAC; iFazio et al. 



20041 ). Observing clusters with a range of times since last disk passage is desirable, but we 
did not use this as a selection criterion because or bital information is only available for a 
fraction of Milky Way globular clusters (e.g., from lOdenkirchen et al.l 119971 ). The clusters 
span a factor of thirty in metallicity (1.5 dex), important for testing the effects of metallicity 



on du st production and retention. All of the clusters are detected with IRAS (IKnapp et al. 



19951 ) , and both NGC 104 (47 Tuc) an d NGC 362 contain stars with known circumstellar 
dust JOridia et al.ll2002l : llta et al.ll2007h . 
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2.1. IRAC observations 

IRAC photometry traces the Rayleigh- Jeans tail of the spectral energy distribution of 
globular cluster stars. Although cold intracluster dust is not expected to radiate at IRAC 
wavelengths, understanding the near-IR SED is very useful for quantifying any IR excess 
from the dust. IRAC observations of the clusters were made over the period 2005 July - 
2006 September as part of Spitzer program identification (PID) 20298 (PI R. Rood), except 
for NGC 5904 which was part of PID 3256 (PI P. Goudfrooij). Dataset identifiers are given 
in Table [2l All but one of the clusters were observed in a single IRAC field of view (5' x 5') 
in all four channels, with 90-252 frames (12 s each) per position (30 frames, 30 s each for 
NGC 5904). NGC 104 was observed in a 5' x 15' map. NGC 104 and NGC 6752 were 
observed in the high-dynamic-range (HDR) mode which includes short (0.6 s) frames taken 
before the 12 s frames. The half-light radius of NGC 104 (190") is slightly too large to 
be fully covered by the PID 20298 observations. Fortunately, observations of this cluster 
covering a wider field of view were available, taken as part of instrument commissioning in 
2003 November (PID 623). These observations were taken in 12-second HDR mode with 3 
frames per position, and cover a 25' x 30' area. 

For all clusters, the Spitzer Science Center (SSC) pipeline version S14 post-BCD mo- 
saics were used for integrated photometry. IRAF0 apphot was used to measure integrated 
magnitudes within the half-light radii given in Table [H For each cluster the background was 
measured in an annulus with inner and outer radii 108 and 117", respectively, and subtracted 
from the integrated magnitudes. For NGC 104 and NGC 6752, the mosaics made from the 
short-exposure frames were used to avoid bright star saturation; the other clusters had no 
saturation problems. The magnitudes measured in the half-light radii were converted to 
total flux densities with multiplication by 2, and corrected with application of the 'extended 
source calibration' given by the SSC. The results are given in Table [21 with electronic links to 
the observational data. The uncertainties in the flux densities due to Poisson noise from the 
cluster and sky background no ise are very small (< 0.5%); uncertainties in the absolute (a 



few percent; iReach et al.ll2005l ) and extended source calibrations (10%) dominate the total 



flux measurements. 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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2.2. MIPS observations 

The MIPS Photometry Astronomical Observing Template (AOT) was used to produce 
5' X 5' images of the cluste rs at 24 and 70 fim. This field size is comparable to the cluster 



sizes (half-light radii from iHarrid [19961, < 3') . Although a somewhat larger field of view 
would have been desirable, this would have substantially increased the observing time. To 
set exposure times, we predicted the 24 a nd 70 fim stellar fiu xes from the clusters' Kg- 



band fiuxes (measured from 2MASS images; iNantais et al.ll2006l ) assuming a Rayleigh- Jeans 
f,y oc z/^ spectrum. The 24 fim observations were required to characterize the cluster SED, 
so high-significance detections of individual stars were not required. The estimated 24 fim 
fiux was used to compute the expected average 24 fj,m surface brightness inside the clusters' 
half-light radii. The resulting values, 0.8-4.2 MJy sr~^, were expected to be detectable in 1 
cycle of the Photometry AOT with 10 second frame times. The intra-cluster dust is expected 
to emit most strongly at 70 fim; however, even at this long wavelength, the integrated fiux 
of the cluster stars is much brighter than the MIPS detection limits. We therefore targeted 
the exposure times to have strong detections of the integrated stellar photospheric fiux; any 
dust signal would appear as an IR excess at 70 fim. The resulting predicted 70 fim fiuxes 
ranged from 6-150 mJy, requiring 1-10 cycles of the MIPS Photometry AOT with 10 second 
frame times. 

The observations were made over the period 2006 June - 2006 November under Spitzer 
PID 30031, a Guaranteed Time program of G. Fazio. The data were processed with the S14 
version of the MIPS pipeline. Post-pipeline processing for the 24 fim BCD images included 
'self-calibration' (dividing by the normalized median image of 'off-source' frames) and mo- 
saicking; the 70 fim BCD images were time- and column-filtere d and then mosaicked. All 



mosaicking used the MOPEX program (IMakovoz fc Khanll2005l ) with the standard param- 



eters recommended in SSC documentation. The 24 fim mosaics have a pixel scale of 2.5", 
and the 70 fim mosaics have a scale of 5.0". Figures [1] and [2] show the cluster mosaic images. 
All of the clusters were detected at high signal-to-noise in the 24 fim images. The same is 
not true at 70 fim, and the major result of this paper is the lack of any strong dust signal 
from any of the eight clusters. 

As for IRAC, fiux densities in the MIPS bands were measured with aperture photometry. 
Inspection of Figure [T] shows that the mid-infrared light for several clusters (NGC 52 72 and 



NGC 6205) appears to have slightly different centers from the optical centers given by [Harris 



(119961 ). We computed new centers by smoothing the 24 fim images and computing the flux 
centroid of the smoothed images in a box of width 4 times the smoothed FWHM, centered 
on the optical center. The mid-IR centers are given in Table [31 Compared to the IRAC 
data, the choice of aperture for the MIPS images is less obvious: the intracluster medium 
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seems more likely to sink to the center of the cluster potential well than to follow the same 
spatial distribution as th e cluster star s . [Ho wever, we note that both the extended emission 
detected in NGC 7078 bvlSover et ali J2006h and the possible intracluster dust in NGC 6024 
reported by iMatsunaga et al.l (120081 ) are offset from the cluster centers.] We decided to use 
two apertures for each cluster: the optical half-light radius as for the IRAC data, and a 
smaller 'core' radius, to maximize signal-to- noise for detection of dust at the cluster center. 
The half-light radii are given in Table [H for these apertures, backg;round annuli with inner 
and outer radii 150 and 175" were used. Following ICarpenter et al.l (120081 ) . we chose a core 
aperture of 16" radius (diameter twice the FWHM of the PSF), with background annuli 
120-140" in radius. To allow for subtraction of stellar contribution to the 70 /im signal 
(see §3.21) . the 24 fim photometry in this aperture was measu red on versions of th e images 
convolved to the 70 fim resolution using a kernel described in iGordon et al.l (120081 ). 



The results of the aperture photometry are listed in Table [3l with electronic links to 
the observational data; no color corrections have been made to these values. The only 
aperture correction made is a mul tiplicative factor of 1.1 3 applied to the large-aperture 
70 /xm photometry. Again following ICarpenter et al.l (120081), uncertainties in the MIPS flux 
density were calculated as 



'7corr^sky''"sky 



sky 

where r^corr accounts for the noise correlation between pixels due to re-sampling during mo- 
saicing (r^corr = 1 for 24 /im and 2 for 70 /xm), ?7sky is discussed below, cisky is the standard 
deviation of the pixels in the background annulus, and iVap and Ns^y are the number of pixels 
in the photometry and background apertures, respectively. (Poisson noise from the clusters 
is negligable compared to the backgro und noise.) After exa mining noise in their 70 fim 
background and photometry apertures, ICarpenter et al.l (120081) concluded that 'excess noise' 
corresponding to r/sky = 1-5 was present in their photometry apertures, a conclusion vali- 
dated by their distribution of signal-to-noise ratios. Although this precise value may not 
apply to our observations, which were made in a slightly different manner, we adopt it here 
to be conservative. The uncertainties in Tab le [3] do not include uncertainties in the absol ute 
cahbration (2% at 24 fim and 5% at 70 fim; iGordon et al.ll2007l : Engelbracht et al.l 120071 ). 



Confusion nois e due to extragalactic sources is not expected to affect these observations: 
Frayer et al.l (120061 ) derived the 70 /xm confusion limit due to extragalactic sources as 0.3 mJy, 
well below the flux regime o f the observa t ions p resented here. Confusion due to Galactic 
cirrus, however, is an issue: IJeong et al.l (120051 ) computed the MIPS 70 /im point source 
detection limits due to Galactic cirrus confusion at the latitudes of our clusters {\b\ > 27) as 
0.1 — 0.2 mJy, again below our detection limits. For larger apertures, however, the confusion 
noise and detection limits scale as d^~°'^'^ where d is the aperture size and a ~ —3.5 is the 
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cirrus power spectrum index. The scaling factors range from about 6 to several hundred 
for the half-light radius apertures used, and hence cirrus confusion will affect our ability to 
detect extended emission for clusters with large angular sizes. 



3. Analysis 



3.1. Cluster spectral energy distributions 

Figure [3] shows the cluster spectral energy distributions (SEDs) from 0.3 to 100 /xm. 
The new IRAC measurements from Table [2] and half-light radius MIPS measurements from 
Table [3] (multiplied by 2 to convert to total fluxes, even though this may not be appropriate 
for non-stellar emission) are shown as filled s ymbols. Also plott ed are extinc t ion-co rrected 
total optical and near-IR magnitudes given by HarrisI (1996) and Cohen et al. (2007). IRAS 



measurements from the Faint Source Catalog (IKnapp et al 



19951 ) are also shown; the 12 and 



25 um measurements h ave been corrected for the cluster extent following the procedure used 
by iKnapp et al.l (119951 ) and flux densities in all IRAS bands have been color-corrected using 
the values for a 5000 K blackbody gi ven in the IRAS Explan atory Supplement. Upper limits 
to the 90 /xm flux densities given by iMatsunaga et al.l (120081 ) are shown for the flve clusters 
in common between the two samples. Overplotted are black-body SEDs derived from flts 
to the flux densities (equally-weighted) at A < 5 /im. All are consistent with the expected 
average temperatures for globular cluster stars (4500-5500 K). 

A flrst check of our measurements is to compare them with previous work. The Spitzer 
and IRAS results are mutually consistent except for NGC 6341, where the 70 /xm data are 
slightly above the IRAS 60 fim upper limits. We suggest that the IRAS limit might be slightly 
too low if there is extended emission from the cluster. For NGC 104, both the 60 and 70 /im 
observations are below the predictions based on the black-body flt, possibly because of the 
large spatial extent o f the cluster. Two of o ur target clusters were observed with the PHOT 
instrument on ISO by lHopwood et al.l (Il999l ): using those authors' method to convert surface 
brightness to flux density, their derived 70 /im flux densities are 130 ± 50 mJy for NGC 104 
and 220 ± 90 mJy for NGC 7252. Our measurements for these clusters are consistent with 
these val ues. The 70 /xm raeasure ments are generally consistent with the 90 /xm upper limits 
given by IMatsunaga et al.l (120081 ) . except for NGC 6341, where the 70 fj,m value appears to 
be too large. The detection for this cluster is at < 3a signiflcance and may be spurious. 

Figure S] compares the predicted (from the black-body flt) and measured flux densities 
(the same values as in the previous flgure) at 24 and 70 fim. We expected that any ICM 
dust would likely be too cool to be detectable at 24 /im; however, circumstellar dust should 
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be warmer and possibly detected at this wavelength. All of the 24 fim flux densities are 
within 20% of the prediction; six clusters are within 10%. This is reasonable accuracy 
considering the difficulties in deriving integrated measurements of these resolved objects. 
The non-detection of 24 fim excesses in integrated photometry may simply reflect the small 
number of mass-losing AGB stars compared to the bulk of the population. Point-spread- 
fitting photometry of the 24 /im images, while beyond the scope of this work, may reveal 
which stars are surrounded by circumstellar dust. Seven of the eight clusters in the sample 
have measurements or upper limits on 70 /im flux that are more than a factor of 2 above the 
black-body prediction. The exception is NGC 104, for which the 70 /xm upper limit is lower 
than the predicted photospheric value. The problem could be background subtraction: the 
half-light radius is almost as large as the 70 fim image and a "background" region is difficult 
to define. 



3.2. Non-stellar fluxes at 70 yum 

Figures [3] and H] clearly demonstrate that only one of the clusters in our sample — 
NGC 6341 — has a significant infrared excess. To quantify this, we computed (limits on) 
the 70 /xm non-stellar emission in both the half-light and core apertures. To compute the 
non-stellar emission in the half-light aperture, we assumed that any photospheric emission at 
70 /xm follows the optical emission, i.e., half of the total emission is contained within R^. As 
an estimate of the stellar emission at 70 /xm we use the predicted value from the black-body 
fit described above. Thus the non-stellar emission at 70 /xm in the half-light radius is 

/70,ns = /70,obs ~ 0.5 X /ycpred- (2) 

The fraction of stellar emission contained within the core aperture is not as straightforward 
to compute; to estimate it, we scale the flux measured at 24 /im in the core aperture (recall 
that the image used to measure this flux was convolved to the 70 fim resolution) by the 
ratios of the predicted total fluxes at 70 and 24 fim. Therefore the non-stellar emission at 
70 fim in the core aperture is: 

/rO.ns = l-76[/70,obs — (/24,obs//24,prcd) X /70,prcd] (3) 

where the factor of 1.76 is the aperture correction for a point source measured in a 16" 
aperture. Table H] gives the derived non-stellar fluxes for each cluster, together with the dust 
mass predictions and observational limits derived in §4.1[ 

Only one cluster shows even tentative evidence for the detection of excess IR emission 
at 70 fim: NGC 6341. Figure [2] shows that the apparent 70 fim emission in this cluster 
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comes from an extended region offset from the cluster core by roughly 43 arcsec, or about 
{2/3)Rh- Within the emission region there are three separate peaks which contain about 
60% of the total flux; it is unclear whether these peaks are separate point sources (none has 
a co-located 24 /^m source) or merely rep resent structure in the emission. In comparison, the 
dust emission detected in NGC 7078 by iBoyer et al.l (120061 ) is much rounder and is located 
closer to th e center of that cluster -Rh/4); however, the possible dust cloud in NGC 5024 
detected by iMatsunaga et al.l (120081 ) is also far from the cluster center, at O.SRh- The IRIS 
60 fim sky maps do show some nebulosity near NGC 6341, so it is possible that the detected 
flux is in fact not associated with the cluster; however NGC 1851, at a similar Galactic 
latitude, also has some nearby nebulosity and does not have a 70 fim detection. Future data 
from the Akari All-Sky Survey should have sufficient spatial resolution to address the issue 
of whether there is diffuse emission near NGC 634 1. The lack of any excess infrared emission 
in the 90 /xm Akari observations of NGC 6341 by IMatsunaga et al.l (120081 ) (see Fig. [3]) casts 
doubt on the reality of our low-significance detection; however, we retain it for analysis in 
the remainder of this paper. 



3.3. Point sources detected at 70 /xm 

While only some of the clusters have spatially diffuse IR emission within Rh, all of 
the 70 /im images except those of NGC 104 and NGC 5904 do contain some faint point 
sources whic h are not concent r ated t oward the cluster centers. [Most of the point sources 



identified by IMatsunaga et al.l (120081 ) are off the fields of view of our images.] What is the 
nature of these sources? The 70 /xm data provide only marginal detections and poor spatial 
localization, but the 24 fim data provide additional clues to the nature of the sources. 

Initially, visual examination was used to generate lists of possible point sources in each 
cluster image. Using an estimated 70 /xm position as a starting point, we performed aperture 
photometry on the 24 and 70 /im images using IRAF/phot, with the photometry routine 
allowed to re-center the sources. Several 70 fim sources were off the edges of the 24 fim 
images, or not visible at the shorter wavelength: the latter group could either be very red 
sources, or noise fluctuations in the 70 /im images. Standard photometric apertures as given 
in the MIPS Data Handbook were used (at 70 /im, aperture radius 18" and background 
annulus 39 — 65"; at 24 /xm, aperture radius 6" and background annulus 20 — 32") and 
appropriate aperture corrections were applied. We did not attempt IRAC photometry for 
the 70 /im sources: the IRAC images of the GCs are signiflcantly more crowded, and we 
were doubtful about the reliability of measurements in these regions. 



The 70 /im sources in the globular cluster flelds could be stars within the cluster, un- 



- 10 - 



related Galactic stars, or background galaxies. To discriminate between these possibili- 
ties, a color-magnitude diagram was made from the MIPS photometry, shown in Figure O 
For comparison, the AGB and post-AGB models (restricted to oxy gen-rich models with 



M < 10 ^Mq, which are most appropriate for comparison with GCs) of lGroenewegenI (120061 ) 



are also shown. The model fluxes have been scaled down by a factor of 5 to reflect a dust-to- 
gas ratio appropriate for GC metallicities. The model fluxes are scaled to a common fiducial 
distance (8.5 kpc) so the same practice was followed for the cluster sources. In addition to 
the evolved star models, we show the (unsealed) flux d ensities and rat i os for 70 fim sources 
in the ELAIS-N2 field of the SWIRE Legacy Survey (ILonsdale et al.ll2003l ): these sources 
are predominantly extragalactic. All but one of the cluster field 70 fim sources are very red, 
with typical values of /70//24 ~ 10, compared to the Rayleigh- Jeans value of ~ 0.1. The 
cluster source colors are more similar to those of the galaxies than to most of the evolved-star 
models, and the stellar models with similar colors are significantly brighter. 

From the colors, we tentatively conclude that the 70 /im sources in the cluster fields are 
likely background galaxies. The hig h Galactic latitude of the clusters (|6| > 26°) supports 
this conclusion, as do the results of Isover et al. ( 2008 ) from their study of 70 /im sources 
near cu Gen and those of Matsunaga et al. (j2008 ). The number density of sources is also 
consistent with this conclusion: there are 23 sou rces in a total of 34 6 arcmin^ for a density 
of 0.066 arcmin"^. The number counts given by lFrayer et al.l (120061 ) for the GOODS-North 
field show that there are 0.07 galaxies arcmin"^ to flux limits o f 11 mJy ( appro priate for 
our average field), comparable to the number in the cluster fields. Ilta et al.l (120071 ) observed 
NGC 104 and 362 at mid-infrared wavelengths from 2.4-24 /im. Those authors found eight 
red sources (F24/F7 > 1) in NGC 362; no su c h sou rces were found in NGC 104. Comparing 
our images of NGC 362 to those of Ilta et al.l (l2007l ). we find that their red sources A, D and 
F are within our field of view. Only source D is detected at 70 /im ; it has -F70/F24 ~ 15, 
consistent with the colors of a galaxy. 



4. Discussion 



4.1. Dust masses and limits 



The non-stellar 70 fj,m fluxes computed in §3.21 allow us to estimate the amount of dust 
present in these clusters. For the clusters wi th negative non-ste llar fluxes, we computed the 
dust mass upper limits as 3crflux- We follow lEvans et al.l (120031 ) in assuming any dust to be 
optically thin, so the dust mass is given by: 



4.79 X 10""/^[mJy] 



© 



D2 



(4) 
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where Ky is the dust absorption coefficient and BuiTd) the Planck function, both in cgs 
units. Our data do not constrain the dust temperature, so we use the value found by 
Bover et all J2006h for NGC 7078 (T^ = 70 K), and the same = 56 cm^ g-\ Different 
assumptions about these t wo quantities can ch ange the derived dust masses by several orders 
of magnitude; for example. lEvans et al.l (120031 ) estimated the dust mass in NGC 7078 by using 
K, for amorphous fayalite at 50 K; at 70 fim this substance has k = 280 cm^ g~\ or fi ve times 
the value used by lBoyer et al.l (120061 ). On the opposite extreme, iKnapp et al.l (119951 ) assume 
silicate gra ins with T = 45 K; wit h cx u, the absorption coefficient at 70 /xm is 36 cm^ g^^. 
Table 5 of iHopwood et al.l (119991 ) shows the effects of different assumptions about the dust 
temperature on derived dust masses: changing from = 70 K to = 40 K, for the 
same composition, increases the inferred by about a factor of ten, while increasing from 
Trf = 70 K to Trf = 100 K decreases by about a factor of 2.5. Uncertainties in the dust 
composition and temperature therefore result in uncertainties of < 2 orders of magnitude in 
derived dust masses. 

Table m gives the dust mass upper limits derived from the 70 /xm measurements in both 
the half-light and core apertures. All values reported are 3a upper limits, except for the 
possible (2.2(t) detection of total dust mass in NGC 6341. This is our only observation 
resulting in even a tentative detection; we show it as a value rather than an upper limit in 
the table and following plots in order to give a visual estimate of the uncertainties. Because 
the cluster emission is generally very faint, we have not attempted to subtract any point 
source emission from the total fluxes; if point sources within the apertures are unrelated to 
the clusters then the true upper limits would be even lower. For all of the clusters in our 
sample, the limits on the total dust masses a re < 4 x 10~'^M, ^ ; the limits on core masses are 



an order of magnitude lower, < 4 x 10 ^Mq. iHopwood et al.l (119991 ) gave upper limits to the 



dust masses for NGC 104 and NGC 5272; our total value for NGC 104 is slightly above their 
e our core valu e is sli ghtly below. Our total value for NGC 5272 is comparable 



limits whi 
to that o f I 

given by iPenny et al. I (Il997h for 80 K dust based on their observations of this cluster at 



Hopwood et al.l (Il999l) for 70 K dust. Our value for NGC 5272 is also below that 



millimeter wavelengths. Our upper limits are above those of iMatsunaga et al. J2008h for the 
five clusters in common, due to the greater sensitivity of their observations. 

The dust mass limits we infer are consistent with measurements and upper limits on 
gas masses. Converting dust masses to gas masses by multiplying by 10^~I^*'/^1, we derive 
gas masses or upper limits for our sample clusters in the ranges 0.002 — 0.0 7Mid (core) and 



0.04- 1.4M^ (total). The only con&med H I detection, O.3M0 in NGC 7078 (Ivan Loon et al. 



20061 ). is nicely within this range. Our inferred limits on the gas mass are consistent with the 
results of prey ious searches for intra-cluster neutral and ionized hydrogen, as su mmarized in 



RobertsI (119881 ). Two clusters have reported detections of ionized gas: NGC 104 (iFreire et al. 
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200l|) and NGC 1851 flGrindlav fc Lilleil 119771 ) . iFreire et al.l (|200l|) derived a total gas mass 



of O.lMfT^ and sta. t ed tha t they expected essentially all the gas in the cluster to be ionized; 



Grindlay fc Lillerl (119771 ) estimated the total mass of ionized hydroge n in NGC 1851 as 



O.O2M0. Our gas mass upper limit for NGC 104 is about th e same as thelFreire et al.l (120011 ) 
detection while our upper limit for NGC 1851 is above the iGrindlay fc Lillen (119771 ) level. 



4.2. Expected dust masses 



To understand the significance of the dust masses and upper limits derived from the 
Spitzer observations, we need theoretical estimates of the expected dust content of the clus- 
ters in the sample. Various estimates of the dust mass within a globular clus ter are possible. 



but lit tle detailed modeling has been done with the exception of the work of lAngeletti et al. 



( 119821 ). Those authors combined dynamical models of star clusters with radiative transfer 
models to predict dust grain temperatures and luminosities for a few individual clusters, given 
the presence of 1O~^M0 of dust. While it would be valuable to repeat such modeling with 
modern radiative transfer codes and understanding of dust properties, for present purposes 
a simp ler estimate of globular cluster dust content will suffice. Following iTayler fc Wood 
( 119751 ) ■ the dust mass may be estimated by: 



Me 



thb 



iVHB(5M) 



100 



(5) 



where thb = 1.8 x 10^ yr is the horizontal branch lifetime, Tc is the time since last crossing of 
the Galactic plane, A^hb is the number of stars on the horizontal branch, {5M) = 0.2Mq is 
the gas mass lost from each star, and [Fe/H] is the cluster metallicity. The last factor scales 
for the change in gas-to-dust ratio with metallici t y. Va lues for Tc are listed in Table [U for 
the two clusters not listed by lOdenkirchen et al.l (119971 ) we assur ne a lower limit of 1 0^ yr . 
To estim ate ./Vhb we use th e 'specific evolutionary fiux' method of iRenzini &: Buzzonil (119861 ) 
(see also lEvans et al.ll2003l ): 

Nub = ^LboiTHB (6) 
assuming B = 2 x 10~^^ st ars yr~^ Lq^ as appropriate for glob ular clusters, and Lboi = 2 x Ly 



(IRenzini fc Buzzonilll986l ) using My and E{B — V) given by iHarrid (119961 ). Fo r one cluster. 



we ca n compare this dust mass estimate to the current dust production rate: lOriglia et al. 



(120071 ) have recently used Spitzer data to identify the mass-losing stars in NGC 104. Summing 
their mass loss rates for the brightest stars (Mboi < —2) and converting to a dust production 
rate with a gas-to-dust ratio of 200, we estimate Md = 1.4 x IO^^Mq yr~^. When multiplied 
by Tc = 6 X 10'' yr, this yields a total dust co n tent o f O.84M0, about twice that predicted 
from Eq. [5l The observations of lOriglia et al.l (120071 ) are likely to be affected by blending; 
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correcting for this would lower the inferred dust mass to a value closer to the prediction. Even 
without such a correction, the comparison gives us confidence that the dust mass predictions 
are of the correct order of magnitude. 

Figure [6] compares the predicted and observed dust masses and upper limits. To aug- 
ment our sample we analyse addi tional clusters obs erved by other gr oups: NGC 635( eI 
faopwood et"al1ll998h . NGC 6656 JPennv et al.lll997h and NGC 7078 JSover et al.ll2006h . 
and the six other clusters observed bv lMatsunaga et all J2008h : NGC 1261, NGC 1904, NGC 
2808, NGC 5024, NGC 5139 NGC 5634 and NGC 6402. Where several values for dust mass 
or upper limits are given in published works, only the one with assumed dust temperature 
closest to 70 K is used. All of our measurements are well below the predicted values; for most 
clusters the disagreement is more than an order of magnitude. The same is true for clusters 
observed by other groups, in that only two c lusters have the inferred dust mass close to the 



predicted value: NGC 6356 (IHopwood et al. 



1998r) fonly fo r specific assumptions about the 



dust distribution), and NGC 5634 (IMatsunaga et al.ll2008l ). one of the lowest-mass, lowest- 
metallicity clusters shown. While assuming a lower dust temperature or higher absorption 
coefficient would increase the dust masses inferred from the observations, it seems unlikely 
that this can explain the three order of magnitude discrepancy for clusters such as NGC 5272. 
The other explanation is that one or more of the assumptions involved in Eq. [5] is incorrect: 
what if globular cluster stars lose much less than O.2M0 on average, the number of horizontal 
branch stars does not scale with Lboi, the dust-to-gas ratio does not scale with metallicity, 
or the dust is removed from clusters on a timescale much shorter than the disk passage time 
Tc? Direct observations of red giant mass loss and c olor-magnitude diagram modeling are 
consi stent with mass losses of approximately O.2M0 ( lOriglia et al.l 120071 : ICaloi fc D'Antona 
20081 ). and the scaling of A^hb with Lboi is well-grounded in stellar evolutionary theory. We 
address the remaining factors below. 

With a small sample of clusters, and an even smaller number of low-significance de- 
tections, deriving conclusions about the effects of GC properties on mass loss is difficult. 
Figure [7] shows the dust masses and limits plotted as a function of [Fe/H], Lboi, Tc, and 
Vpf^r^n- (In this and the following figure, we plot the more stringent dust mass limits of 
Matsunaga et al.l (120081 ) for the five clusters in common with our sample.) There are no 
clear correlations: while the clusters with large inferred dust masses tend to have large Tc, 
consistent with the expectation that intracluster dust builds up between cluster passage 
through the disk, other clusters with similar values of Tc have much smaller inferred masses. 



^For this object we plot the more conservative 'dust in beam' mass given by iHopwood et al.l (|l998f ): 
Matsunaga et al.l (j2008l ) use the larger 'dust in core' value. Regardless of this choice, this cluster still has 
the largest inferred mass. 
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This agrees with the finding of iBoyer et aLl (120061 ) that the ICM dust in NGC 707 8 must 
have been collecting for a time t <^ r^, and with the computation by iMatsunaga et al.l (120081 ) 
that cluster dust lifetimes -C Tc. The clusters with dust detections do not have the highest 
central escape velocities, suggesting that possession of a deep potential well is not in itself 
enough to guarantee the presence of an ICM. The clusters cover a relatively small range in 
bolometric luminosity; this parameter does not seem to have a strong effect on inferred dust 



mass. 



How does metallicity relate to intracluster dust? The two clusters with the largest 
inferred dust masses (NGC 7078 and NGC 6356) have the lowest and highest metallicities 
in the sample. NGC 5024, with a tentative dust detection, also has a low metallicity. The 
globular cluster data do not help to clarify the controversy over the connection between 
metallicity and dust production: while some studies show that met allicity and mass loss are 



not correlated in carbon- rich AGB stars in the Magellanic Clouds (ILagadec fc Ziilstrall2008 



Sloan et al.l 2008), others fin d that metal-poor carbon stars, like oxygen-rich giants, form less 
dust (Ivan Loon et al.l 120081 ). The clusters considered here are almost all more metal-poor 
than the Magellanic Clouds, so their AGB stars are expected to be oxygen, rather than 
carbon-rich; however the effects (if any) of metallicity on dust mass loss are not evident from 
the present data. 



4.3. Dust removal mechanisms 



Observations of mass-losing stars within GCs show that mass loss is occurring at roughly 
the expected rate, so the non-detection of the ICM means that it must somehow be removed. 
One possible mechanism is ram-pressure stripping of the cluster ICM by the Galactic halo. 
Support for this mechanis m is provided by detections of diffuse X-ray emissio n near several 



Galac tic globular clusters (lOkada et al.ll2007l : iKrockenberger &: Grindlaylll995l ): lOkada et al. 



(120071 ) show that the temperature and luminosity of the diffuse emission is consistent with it 
arising from stripped intracluster gas, and that the energetics require nearly all of the cluster 
gas to be removed. Those authors also state, to show diffuse X-ray emission, a cluster should 
have a velocity relative to the Galactic halo > 150 km s~^. We corn p uted Vh for our sampl e 
cl usters from the space velocit ies given in lOdenkirchen et al.l (119971 ) , iDinescu et al.l (119991 ) , 



or 



Casetti-Dinescu et al.l (120071 ) assuming a non-rotating halo; these are plotted against dust 



masses and limits in the left panel of Figure [HI The high velocities of both NGC 7078 and 
NGC 5024 appear to be problematic for halo stripping; a measurement of the space velocity 
of NGC 6356 might help to clarify the situation. There is certainly no evidence that clusters 
with low Vh have more dust; however, the correlation of cluster kinematics with metallicity 
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combined with any dependence on metallicity of the dust mass hmits may also confuse the 
issue. 

Another po ssible ICM rem o val me chanism is kinetic energy injected by stellar collisions, 
as proposed by lUmbreit et al.l (120081 ) for M15. Those authors use an N-body model to 
estimate the time between collisions in the cluster core; their result is of the same order of 



magnitude as a less-sophisticated estimate based on Eq. 8-125 of iBinney &: Tremaind (119871 ): 



^coU 



[7.9 X lO'Vr] 



plRl 



(7) 



where N is the number of stars within the cluster core (of radius Rc pc), Uc is the cen- 
tral velocity dispersion in km s~^, and po is the central mass density in M0 pc~'^. Using 
the values for -R^, pn, and <Jr = ('2 7T)~^^'^rhon^^ fn from from King-model fits given in 



McLaughlin &: van der Marell (120051 ) . we computed tco\\/N for each cluster in the sample; it 
is plotted against in the right panel of Fig. [SI Again, there is no clear correlation of the 
time between collisions and dust masses; while both NGC 6356 and NGC 5024 have long 
tcoii, so do many other clusters. NGC 7078 is core-collapsed and should have a short central 
collision time. We conclude that neither ram pressure stripping and stellar collisions as the 
cause of ICM removal can be ruled out by the current data. 



5. Conclusions 



The intracluster medium of globular clusters remains elusive. Our sample comprised 
eight Galactic globular clusters thought to have the best chance for retaining an intracluster 
medium. No significant extended emission is detected from seven of the eight target clusters: 
any such emission may well be hidden underneath confusion noise from Galactic cirrus. One 
cluster, NGC 6341, has tentative evidence for excess IR emission from cool dust, although the 
lack of a 90 pm excess casts doubt on the detection. All of our inferred dust masses and upper 
limits are < 4 x IO^^Mq, compared to predicted dust masses ranging from 10^°'^ — 10^^'^ Mq 
for the same clusters . Searches for intracluster gas have been similarly unsuccessful (e.g., 
van Loon et al.ll2006l ). The lack of intracluster dust and gas implies that either (1) evolved 



stars in globular clusters make less dust than predicted, or (2) the ICM escapes from the 
clusters more quickly than expected. The first possibility seems unlikely because of the 
evidence to the contrary in both NGC 7078 and NGC 5139, where t wo of the three dustiest, 
most mass-losing stars have [Fe/H] < —2.0 (IBoyer et al.ll2006l . l2008l ). 



It is important to point out that ICM escape is not just a trivial occurrence in the life 
of a globular cluster: substantial mass loss over a GC's history could affect its dynamical 
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evolution and perhaps even its survival as a bound system (Ivan Loon &: McDonaldl l2006l ) . 
Several mechanisms for ICM escape have been discussed above, with energy injected by 
stellar collisions perhaps being favored over ram pressure from Galactic halo gas. Stellar 
winds provide another possible mechanism; gia nt stars with wind velocities above th e escape 
velocity of a typical cluster have been observed ( Smith et all2004 : Dupree et al. 1992 ). While 
there is observation al evidence that metallicity does not strongly a ffect the mass loss rate 
of individual stars (IMcDonald &: van LoonI 120071: ISloan et al.ll2008l ). stellar wind velocities 



do appear to scale with metallicity ( Marshall et al.ll2004l ). Faster winds in AGB stars in 
metal-rich clusters could therefore increase the rate at which the ICM is driven out of the 
cluster. The existing small sample of clusters does not permit us to distinguish between the 
possibilities for ICM removal. 

Decades of searching have yielded few detections of the elusive globular cluster ICM. A 
possible strategy for future observations hoping to detect the ICM is to study the removal 
mechanisms in detail and identify those clusters most likely to retain some ICM. Observa- 
tions of distant clusters which have had few interactions with the Milky Way could place 
limits on the importance of halo stripping, although the larger distance of such clusters re- 
duces the effective sensitivity. While infrared facilities available in the near future will not 
greatly improve on the sensitivity of Spitzer, observing a larger sample of clusters would 
allow a more thorough exploration of the dust content as a function of cluster properties. 
Observations with future radio facilities such as the Square Kilometer Array will place more 
stringent limits on the gas content. Together with more sophisticated modeling efforts and 
detailed examination of mass loss in individual cluster stars, a multi-wavelength observa- 
tional approach provides the best chance for finding and explaining the elusive intracluster 
medium in globular clusters. 
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Fig. 1. — MIPS 24 /iin images of Galactic globular clusters. Each image is 5.2' across, with 
north up and east to the left. The black or white circles show the half-light radii. 




Fig. 2. — MIPS 70 yum images of Galactic globular clusters. Each image is 5.2' across, with 
north up and east to the left. The black circles show the half-light radii (values given in 
Table ^ 
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Fig. 3. — Spectral energy distributions of Galactic globular clusters, from 0.3 to 100 /xm. 
Filled symbols are IRAC and MIPS measurements from this work. O pen syrnbols a nd grey 
upper limits are pu blished measuremeri ts: optical magnitudes from iHarrid (1l996l ). near- 
IR magnitud e s from ICohen et al.l (120071 ). IRAS F aint Source Catalog me asurements from 
Knapp et al.l (Il995l ). and 90 /im upper limits from iMatsunaga et al.l (120081 ). 
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Fig. 4. — Comparison of predicted and measured flux densities for globular clusters at 24 /xm 
(squares) and 70 /xm (stars and upper limits). 
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Fig. 5. — MIPS color-magnitude plot for sources detected at 70 /xm (squares ), AGB star 
(open circles) and post- AGB star (filled triangles) models from iGroenewegenI (120061 ). and 
SWIRE ELAIS-N2 sources (crosses). The AGB models and cluster source fluxes have been 
scaled to a distance of 8.5 kpc, assuming the sources to be at the distance of the cluster. 
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Fig. 6. — Cluster dust masses predicted with Eq. [5] versus derived masses and upper limits 
from the MIPS observations for both total dust mass (filled black symbols) and core dust 
mass (filled grey symbols). Previously published dust limits are open symbols, not labeled 
with cluster NGC number. Dashed line is the line of equality. 
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Fig. 7. — Derived cluster dust masses and upper limits versus time since last Galactic 
disk passage (upper left), central escape velocity (upper right), metallicity (lower left), and 
bolometric luminosity (lower right). Filled symbols: results from this work (total masses), 
open symbols: previously published work. An uncertaint y of 20% has beeri assig ned to Md 
for NGC 5024 as no value was given in the original work (IMatsunaga et al.ll2008l ). 
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Fig. 8. — Derived cluster dust masses and upper limits versus velocity relative to the halo 
(left) and estimated time between stellar collisions in the core (right). Symbols as in Fig. [71 In 
the left panel, clusters NGC 1261, 5646, 6356, and 6402 are not shown, as orbital parameters 
were not available. In the right panel, b oth NGC 6752 and NGC 7078 are plo tted at tcou/^ = 
10^ yr; these clusters are core-collapsed (IMcLaughlin fc van der Marelll2005l ) and should have 
short collision times. 
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Table 1. Globular cluster targets 



Target 


Distance^ 














kpc 




km s"-'- 


arcsec 


Lq 


yr 


NGC 104 (47 Tuc) 


4.5 


-0.76 


65.9 


190.1 


1.12 X 10"^ 


6 X 10^ 


NGC 362 


12.1 


-1.22 


46.7 


49.1 


4.54 X 10^ 


3 X 10'' 


NGC 1851 


8.5 


-1.16 


48.7 


30.4 


3.88 X 10^ 




NGC 5272 (MS) 


10.4 


-1.57 


27.6 


138.4 


6.56 X 10^5 


2.1 X 10** 


NGC 5904 (M5) 


7.5 


-1.27 


37.2 


105.9 


6.21 X 10^ 


2.1 X 10** 


NGC 6205 {M13) 


7.7 


-1.54 


32.1 


101.4 


5.46 X 10^ 


7 X 10' 


NGC 6341 {M92) 


8.2 


-2.28 


35.9 


61.4 


3.44 X 10^5 


6 X 10' 


NGC 6752 


4.0 


-1.56 




114.8 


2.36 X 10^ 





''Distance, metallicity and bolometric luminosity from iHarrij l|l996h . 

^Proj ected half-light radius and centr al escape velocity tiesc.o from Ki ng-model fits 
given in lMcLauehlin fc van der Marell ||2005| ') (except i?,,, for NGC 6752, from llYaeer" et al.l 
1 19951 '). 

"^Time since last passage through Galactic disk Tc, from lOdenkirchen et al] lll997l ). 



Table 2. IRAC integrated photometry 



Target 


/,(3.6) 


^(4.5) 


/,{5.8) 


/.{8.0) 


Datasct ID 




Jy 


Jy 


Jy 


Jy 




NGC 104 


107 


67.9 


46.4 


24.7 


ADS /Sa. Spitzor#0014502656, ADS/Sa. Spitzer#0007860992 


NGC 362 


7.10 


4.82 


3.44 


1.73 


ADS/Sa. Spitzer#0014503168 


NGC 1851 


3.57 


2.31 


1.64 


0.82 


ADS/Sa.Spitzer#0014503424 


NGC 5272 


8.67 


5.68 


3.63 


1.85 


ADS/Sa.Spitzcr#0014504192 


NGC 5904 


9.25 


6.56 


4.81 


2.42 


ADS/Sa.Spitzer#0011586304 


NGC 6205 


9.66 


6.43 


4.35 


2.17 


ADS/Sa. Spitzer#0014504704 


NGC 6341 


5.05 


3.38 


2.31 


1.16 


ADS/Sa.Spitzer#14504960 


NGC 6752 


13.25 


8.41 


6.32 


3.42 


ADS/Sa.Spitzer#14506240 
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Table 3. MIPS cluster photometry 



half-light core Dataset IDs 

Target Center /^(24) /^(70) /^(24) /^(70) 

[J2000.0] mjy mjy mjy mjy 



NGC 104 


00:24:07.4 -72:04:48 


1980 ± 18 


124 + 186 


91 + 2 


NGC 362 


01:03:14.8 - 70:50:51 


138 + 5 


42 + 22 


29 + 1.2 


NGC 1851 


05:14:06.8 -40:02:49 


63 + 3 


-7 + 11 


20 + 1 


NGC 5272 


13:42:11.6 +28:22:47 


161 + 6 


1 + 87 


10 + 0.7 


NGC 5904 


15:18:33.0 +02:04:50 


172 + 6 


-34 + 67 


14 + 0.8 


NGC 6205 


16:41:39.9 +36:27:28 


153 + 5 


67 + 62 


9 + 0.6 


NGC 6341 


17:17:07.8 +13:08:11 


88 + 4 


68 ± 25 


13 ± 0.8 


NGC (iTr,2 


i<):ll):r,2.0 -o9:59:()2 


218 ± (i 


-.') ± 129 


13 ± 0..S 



2 + 9 ADS /Sa.Spitzer#0017315328, ADS /Sa.Spitzcr#001731353f 

5 + 5 ADS/Sa.Spitzer#0017316352, ADS/Sa.Spitzer#001731456( 

-1 + 4 ADS/Sa.Spitzer#0017316096, ADS/Sa.Spitzer#001731404J 

-2 + 5 ADS/Sa.Spltzer#0017316608, ADS/Sa.Spltzer#00173148ie 

-9 + 5 ADS/Sa.Spitzer#0017313024, ADS/Sa.Spitzer#001731276J 

-1 + 6 ADS/Sa.Spitzcr#0017315584, ADS/Sa.Spitzcr#001731379^ 

-1 ± 5 ADS/Sa.Spitzcr#001731."5840. ADS/Sa.Spitzcr#001731'i30- 

-1 + 9 Al)S/Sa.Si)i(z('r#(K)17315072. ADS /Sa. Spit zc-r#00173i328( 



Note. — All uncertainties are la. 



Table 4. Non-stellar 70 /xm fluxes and inferred dust masses 



Target 


/to, tot 
mJy 


/yOjCorc 

mJy 


M© 




•M"total 
M© 




.^^core 

M© 




NGC 104 


-75 + 186 


-15 + 9 


4.7 X 10" 


-1 


< 1.5 X 10- 


-4 


< 1.2 X 10" 


-5 


NGC 362 


26 + 22 


3 + 5 


3.3 X 10" 


-2 


< 6.3 X 10" 


-5 


< 2.5 X 10" 


-6 


NGC 1851 


-14 + 11 


-6 + 4 


> 1.1 X 10" 


-3 


< 6.4 X 10" 


-5 


< 4 X 10" 


-6 


NGC 5272 


-15 + 87 


-6 + 5 


1.5 X 10" 


-1 


< 3.7 X 10- 


-4 


< 3.5 X 10- 


-5 


NGC 5904 


-54 + 67 


-18 + 6 


2.8 X 10" 


-1 


< 1.0 X 10- 


-4 


< 1.5 X 10" 


-5 


NGC 6205 


48 + 62 


-4 + 6 


4.4 X 10" 


-2 


< 1.5 X 10- 


-4 


< 2.3 X 10" 




NGC 6341 


59 + 25 


-4 + 5 


4.3 X 10" 


-S 


(5.2 + 2.2) X 10- 


-5 


< 2.2 X 10" 


-5 


NGC 6752 


-28 + 129 


-10 + 9 


> 2.6 X 10" 


-4 


< 8.2 X 10- 


-5 


< 9.5 X 10- 


-6 



